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ABSTRACT 

The X-ray properties of groups and clusters of galaxies obey scaling relations that 
provide insight into the physics of their formation and evolution. In this paper, we 
constrain gas and dark matter parameters of these systems, by comparing the observed 
relations to theoretical expectations, obtained assuming that the gas is in hydrostatic 
equilibrium with the dark matter and follows a polytropic relation. In this exercise, we 
vary four parameters: the gas polytropic index F, its temperature at large radii, the 
■ dark matter logarithmic slope at large radii C, and its concentration. When comparing 

the model to the observed mass-temperature relation of local high-mass systems, we 
find our results to be independent of both the gas temperature at large radii and 
of the dark matter concentration. We thus obtain constraints on T, by fixing the 
dark matter profile, and on ^, by fixing the gas profile. For a Navarro-Frenk- White 
dark matter profile, we find that F must lie between 6/5 and 13/10. This value is 
consistent with numerical simulations and observations of individual clusters. Tak- 
ing 6/5 F < 13/10 allows the dark matter profile to be slightly steeper than the 
Navarro-Frenk- White profile at large radii. Upon including local low-mass systems, we 
obtain constraints on the mass-dependence of F and on the value of the gas temper- 
ature at large radii. Interestingly, by fixing F = 6/5 and C, = —3, we reproduce the 
observed steepening/breaking of the mass-temperature relation at low masses if the 
temperature of the intercluster medium is between 10^ K and 10*" K, consistent with 
numerical simulations and observations of the warm-hot intergalactic medium. When 
extrapolated to high redshift, the model with a constant F reproduces the expected 
self-similar behaviour. Given our formulation, we can also naturally account for the 
observed, non-self-similar relations provided by some high-redshift clusters, as they 
simply provide constraints on the evolution of F. In addition, comparing our model to 
the observed luminosity-temperature relation, we are able to discriminate between dif- 
ferent mass-concentration relations and find that a weak dependence of concentration 
on mass is currently preferred by data. In summary, this simple theoretical model can 
still account for much of the complexity of recent, improved X-ray scaling relations, 
provided that we allow for a mild dependence of the polytropic index on mass or for 
a gas temperature at large radii consistent with intercluster values. 
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1 INTRODUCTION 

Clusters of galaxies are important probes of cosmology, 
as they complement other observations, including Type la 
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supernovffi (e.g. Riess et al. 1998; Perlmutter et al. 1999) 
and the cosmic microwave background (e.g. Komatsu et al. 
2009), in building the currently favoured "concordance cos- 
mology" model. As the largest gravitationally bound sys- 
tems in the Universe, galaxy clusters trace the growth of 
structure, which in turn strongly depends on the cosmolog- 
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ical parameters (for recent reviews on clusters in a cosmo- 
logical context, see e.g. Voit 2005; Allen, Evrard & Mantz 
2011). However, galaxy clusters can provide useful cosmolog- 
ical information only if their mass, composed predominantly 
of dark matter (DM), can be measured accurately enough. 
Given the elusive nature of DM, this can be achieved only 
through indirect methods, e.g. by measuring orbital veloc- 
ities of galaxies (e.g. Zwicky 1933, 1937), by gravitational 
lensing techniques (e.g. Zwicky 1937), or through studying 
the intracluster medium (see reviews by e.g. Sarazin 1986; 
Bohringer & Werner 2010). 

The intracluster medium is the hot, X-ray emitting gas 
that to first order is in thermal equilibrium with the clus- 
ter gravitational potential. One of the most commonly used 
models for galaxy clusters utilises a King (1962) profile to 
describe the total gravitational potential with the isothermal 
gas in thermal equilibrium. Solving the hydrostatic equilib- 
rium (HE) equation, the resulting profile, referred to as the 
/9-model (Cavaliere & Fusco-Femiano 1976, 1978), is widely 
utilised to describe radial profiles of clusters. However, de- 
tailed X-ray observations of individual clusters have shown 
clear radial gradients in their temperature profiles, hence 
calling into question the assumption of isothermality. 

One of the simplest extensions of the isothermal model 
is to assume a poly tropic relation, wherein gas density and 
pressure (or temperature) are related by a simple power law. 
This description naturally provides radial temperature gra- 
dients. A polytropic relation is indeed expected in ideal adi- 
abatic processes, with the polytropic index in those cases 
being equal to the ratio of specific heats at constant pres- 
sure and constant volume. This value of the exponent, how- 
ever, does not appear to be consistent with observations (e.g. 
Bode, Ostriker & Vikhlinin 2009 and references therein), 
suggesting that the polytropic relation is not an equation of 
state. Although the origin of the observed polytropic rela- 
tion is still a puzzle (but see Bertschinger 1985 for a possible 
explanation), HE-polytropic solutions have been considered 
since the early 1970s (e.g. Lea 1975). Hybrid descriptions 
have also been considered, where the polytropic relation and 
the /3-model description for the gas density are employed si- 
multaneously (Markevitch et al. 1998; Ettori 2000; Sander- 
son et al. 2003; Ascasibar et al. 2003). 

Numerical prescriptions of DM have advanced to the 
point where we are fairly certain about the profile of the DM 
potential (e.g. Navarro, Frenk & White 1996: NEW). Suto, 
Sasaki & Makino (1998) first derived the HE-polytropic so- 
lution in the case of a gravitational potential described by 
an NFW profile. The concept is simple: once the DM has 
been completely described (in the case of NFW, by only 
two parameters, e.g. virial mass and concentration), the gas 
is analytically described with just three free parameters: the 
polytropic index and the values of two gas quantities, cho- 
sen amongst e.g. density, pressure and temperature, at some 
given radius (usually zero or infinity). Extracting these three 
gas parameters and checking results with observations is not 
entirely straightforward. Radial profiles of individual clus- 
ters (e.g. Pratt et al. 2007; Arnaud et al. 2010) can provide 
some hints, but resolution limits and scatter prevent us from 
obtaining accurate constraints on the three unknown scalars 
introduced above. Moreover, any uncertainty in the DM pro- 
file makes the constraining more difficult. More global scal- 
ing relations between cluster mass, temperature and lumi- 



nosity, on the other hand, can be very helpful. If one has a 
model for the temperature profile, then one can also calcu- 
late the average temperature and obtain a theoretical mass- 
temperature (M-T) relation that can be easily compared to 
observations (e.g. Vikhlinin et al. 2006; Vikhlinin et al. 2009: 
V09; Sun et al. 2009: S09). This is what we largely pursue in 
this paper. Since these average quantities are usually com- 
puted after excising the central region of the cluster, any 
uncertainty in the properties of the inner parts of the DM 
halo becomes unimportant. On the other hand, any change 
in the outer regions of the cluster can produce significant 
effects on the average quantities. For this reason, we do not 
limit ourselves to the NFW profile, but use its generalised 
version by Bulbul et al. (2010: BIO), in which the DM log- 
arithmic slope at large radii is not fixed. For completeness, 
given the uncertainty in the mass-concentration relation, we 
also let the concentration vary, thus having two DM free 
parameters. 

The choice of the three primary gas quantities and the 
method for constraining them vary between authors. Ko- 
matsu & Seljak (2001: KSOl), in their description of a uni- 
versal gas profile, assume that the gas density profile is equal 
to that of DM in the outer regions of the DM halo (assumed 
to be described by an NFW profile), obtaining a constraint 
on both the polytropic index and the central gas temper- 
ature. Their assumption is reasonable and very powerful, 
as it gives constraints on the polytropic index itself, which 
has not been done in other studies. However, it predicts an 
M-T relation which is not in good agreement with recent 
X-ray observations. Also, the gas profile is not completely 
described, as it needs a second constraint to fix the central 
density. Komatsu & Seljak (2002: KS02) provide such a con- 
straint by requiring that the gas density in the outer regions 
of the DM halo is equal to the DM density multiplied by the 
cosmic baryon fraction. 

Ascasibar et al. (2003), investigating the differences be- 
tween polytropic relations and /^-models, fix instead the 
polytropic index as an output from their numerical simula- 
tions and obtain constraints on the central gas temperature 
and density, by imposing the gas density to be zero at in- 
finity and the baryon fraction to never exceed the cosmic 
value. Ascasibar et al. (2006: A06), studying the origin of 
self-similar scaling relations, use a similar method, with the 
slight difference that the baryon fraction at three times the 
DM scale radius must be equal to the cosmic value. 

More recently, a series of papers (Ostriker, Bode & 
Babul 2005; Bode et al. 2007; Bode et al. 2009) have ex- 
plored in detail models with polytropic gas in HE with a 
DM potential, adding additional physics such as star forma- 
tion and feedback, in the form of simple recipes. Similar to 
the approach of Ascasibar et al. (2003) and A06, Ostriker et 
al. (2005) take the polytropic index as a given (from numeri- 
cal simulations and observations) and additionally constrain 
the central gas pressure and density by imposing conserva- 
tion of energy and by matching the external surface pressure 
to the momentum flux from the infalling gas at the virial ra- 
dius. The second constraint arises from a sharp truncation 
in the DM proflle at the virial radius, adopted by Ostriker et 
al. (2005), but not contemplated in KSOl, KS02, Ascasibar 
et al. (2003), A06, nor in this paper. By adding terms to 
the energy conservation equation, they are able to discrim- 
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inate between simple feedback recipes when matching the 
resulting models to observations. 

In this paper, we assume that an ideal polytropic gas is 
in HE with the DM in groups and clusters of galaxies and 
explore the consequences of the polytropic solution, consid- 
ering different ways to specify the three primary gas quanti- 
ties and the two free DM parameters. This is done assuming 
no prior knowledge of either the polytropic index or the 
boundary conditions, and only partial knowledge of the DM 
profile. Thus, our approach differs from previous work in the 
following ways, (i) We do not assume the DM to be NFW- 
like at large radii; instead, we adopt a generalised version 
of the NFW profile (BIO), and are able to discriminate be- 
tween different values of the logarithmic slope of the DM at 
large radii, (ii) We do not fix the value of the polytropic in- 
dex from simulations (Ascasibar et al. 2003; A06; Ostriker et 
al. 2005; Bode et al. 2007, 2009) nor constrain it with some 
assumption on the gas profile at large radii (KSOl; KS02); 
instead, we leave the polytropic index as a free parameter 
and, by comparing our model to observations, constrain it 
as a function of mass and redshift. (iii) We do not force the 
gas radial profile to behave in any way, either by imposing 
its profile to be equal to that of DM at large radii (KSOl; 
KS02) or by requiring it to vanish at infinity (Ascasibar 
et al. 2003; A06); instead, we simply require the temper- 
ature of the gas to never be negative. By leaving the gas 
temperature at large radii as a free parameter, we are able 
to constrain the value of the temperature of the intracluster 
medium (or of the intercluster medium, depending on the as- 
sumptions we make) by matching our theoretically derived 
M-T relation to observations in the low-mass regime (i.e. 
groups), (iv) We do not fix the concentration from simula- 
tions (KSOl; KS02; Ascasibar et al. 2003; A06; Ostriker et 
al. 2005; Bode et al. 2007, 2009), but instead investigate the 
effects of different mass-concentration relations; we discrim- 
inate between different recipes by matching our model to 
recently observed luminosity-temperature (L-T) relations 
(Maughan et al. 2011: Mil). For this last part of the study, a 
second constraint was necessary, and we imposed the baryon 
fraction within the virial radius to be equal to the cosmic 
value. 

We believe the flexibility of our model to be its main 
strength. More general assumptions than in previous studies 
help us obtain more constraints: on the polytropic index; on 
the DM logarithmic slope at large radii; on the gas temper- 
ature at large radii; on concentration. Also, we do extend 
the model to higher redshift and obtain constraints on pa- 
rameter evolution. 

Throughout the paper, we use the following cosmologi- 
cal parameters (WMAP5; Komatsu et al. 2009): Qo = 0.258, 
fib = 0.0441, = 0.742, Ho = WOh km s"^ Mpc"^ = 
70/170 km s"^ Mpc"\ h = 0.719 (i.e. /170 = 1.027). The val- 
ues and uncertainties of other parameters (e.g. erg) are not 
relevant in our study. 

The outline of the paper is as follows. In Section [2l we 
formulate the problem and study in detail the consequences 
of assuming hydrostatic equilibrium and a polytropic rela- 
tion. In Section[3l we derive galaxy groups and clusters scal- 
ing relations and compare them to recent observations to 
obtain constraints on gas and DM parameters. We conclude 
in Section [4] 



2 POLYTROPIC GAS IN GALAXY GROUPS 
AND CLUSTERS 

2.1 Formulation 

In this section, we model groups and clusters of galaxies as 
isolated, spherically symmetric, composite systems of stars, 
DM and ideal gas. We describe the DM with a generalised 
NFW density profile (BIO), 

p-DM(r) = W 

(r/r-s)(l +r/rs)P 

where 1 < /? < 3, rs is the scale radius, Sc is the charac- 
teristic (dimensionless) density, Pc{z) = 3Hq E'^ (z) / (SnG) is 
the critical density of the Universe at redshift z, E{z) = 
H{z)/Ho = [no{l + zf + (1 - Qo - Qa){1 + zf + J^a]'/' 
and G is the gravitational constanlQ. Note that the loga- 
rithmic slope at small radii is the same as for the standard 
NFW profile, which is a particular case of this profile when 
(3 = 2. We further assume the gas to follow a polytropic 
relation, with the pressur^l P, density p and temperature 
T related as follows: P{r) = Po[p{r)/po]^ or, equivalently, 
p{r) = po[T{r)/To]^^'-^'^\ where Po, po and To are the 
central values of pressure, density and temperature, respec- 
tivelv, and ci(log P) /d(log p) = F > 1 is the polytropic in- 
dejQ, assumed to be the same at all radii (i.e. the polytropic 
relation is complete). 

We assume the gas to be in global, thermally-supported 
HE with the total gravitational potential 4>, that is, to obey 
the HE equation, VP = —pVcj), everywhere. We also neglect 
the effects of gas self-gravity, as it has been demonstrated 
(e.g. Suto et al. 1998) that the contribution of the gas to 
the total gravitational potential does not lead to signifi- 
cantly different results, and additionally ignore the stellar 
gravitational potential, as it has been shown (see Capelo, 
Natarajan & Coppi 2010 for the low-mass case) that gas 
profiles are affected by stars only in the very inner regions, 
which are not important for the results of this worljj. The 
total gravitational potential is thus dominated by that of 
DM, i.e. (j) = ^DM, which is given by 

0DM(r)= <^ ^V-2(r/r.)(l + r/rs)f-2 " ^ ^' (2) 
Uoln[l-h(r/r,)]/(r/r,) if /? = 2, 

where 4>o = 0dm (0) = — 47rG(5cPc'"s/(/3 — 1) and we have 
used Equation ^ to solve Poisson's equation, V^0dm ~ 
AnGpDM, after imposing the Dirichlet boundary condition 
0dm(+oo) = 0. 

^ This profile and /3 should not be confused with the /3-model of 
Cavaliere & Fusco-Femiano (1976, 1978). 

^ Unless otherwise stated, pressure is thermal gas pressure, den- 
sity is gas mass per unit volume, and temperature is gas temper- 
ature. 

^ In other notations, the polytropic index n is related to F via 
n=l/(r-l). 

* To be precise, the stellar gravitational potential has a strong 
effect on the central gas quantities (e.g. Capelo et al. 2010), but 
not on the luminosity and average temperature, when calculated 
with an excision of the central region (see Section [Sj. For the 
same reason, we do not investigate DM profiles with different 
logarithmic slopes at small radii (e.g. Moore et al. 1998). 
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The solution to the HE equation is then given by (Suto 
et al. 1998) 



p{r) = po 



r/ra 



(3) 



for (3 = 2 and by 



P{r) = po 



r - 1 « 
1 A„ 



(4) 



(5) 



for /3 7^ 2, where the dimensionless gas parameter 

_ po _ ^nGScPcrsfJ-mp 

ks is the Boltzmann constant, rrip is the proton mass and 
p = 0.62 is the mean molecular weight of the gas (e.g. Par- 
rish et al. 2011), assumed to be completely ionized and to 
have mean metallicity [Fe/H] = —0.5. 

In both cases, we used the boundary condition p(0) = 
po, where po can in principle have any non-negative value 
(although we will see later in this section how to con- 
strain it). BIO imposed a more stringent boundary condi- 
tion, T(-l-oo) = Too ~ 0, and obtained a solution which is 
a particular case of Equation (U) when (F — l)Agas/F — 1. 
Our solution is a more general case that allows for any non- 
negative value of Too , which can be then fixed depending on 
the physical situation (see Section [3. 1.21 for more details). 

We now specify the problem further and consider a 
system with a given virial mass Mvir at a given redshift 
z, assuming that the virialisation was reached right before 
the time at which we observe it (recent-formation approx- 
imation). In this paper, unless otherwise stated, we define 
virial mass as the total (stars, DM and gas) mass within 
the virial radius, defined below. The DM density profile 
given in Equation ([l]) can then be written, via the equal- 
ity 5c = /DMAvirCvir/[3/(Cvir)], aS 



PDM(r) = 



A/DM(rvir) 



(6) 



where the virial radius rvir = [3Mvir/(47rAvir/5c)]^'''' and the 
concentration Cvir ~ '"vir/j's are both assumed to be inde- 
pendent of /3; Avir is the virial overdensity, described by 
Avir(z) ^ IStt^ + 82[Qm{z) - 1] - 39[nn,{z) - if (Bryau 
& Norman 1998); the matter fraction r2m(z) = S^o(l + 
zf/E'^{z); /dm = MDM(rvir)/Mvir = 1 - .fb is the virial 
(that is, within the virial radius) DM mass fraction; /b = 
Qh/^o is the virial baryon mass fraction, set equal to the 
cosmic value, assumed constant with redshift and mass; the 
function f{x) = ln(l + x) — x/(l + x), when /3 = 2, or 
fix) = [l/(/3 - 1) + [-X + 1/(1 - /3)]/(l + xf-']/{p - 2), 
when /3 7^ 2. We note that several authors use different 
definitions for mass and radius. In general, one can define 
such that Ma^ = 47rrA^ AcPc/3, where Ac is the over- 
density ratio. This ratio can be either a fixed number (e.g. 
200, 500, etc.) or can be replaced by the just defined Avir, 
so that TAc is the virial radius rvir and Ma^ is the virial 



mass Mvir. See Appendix \K\ for the calculation of the re- 
lations between different overdensity masses and radii, for 
both the NFW case and for the more general BIO case. 
Moreover, some authors use equivalent but different nota- 
tions for the virial overdensity: Afvir ~ 47rrvir AvirJ^m(z)pc/3 
and Avir(2) ^ [18n^+S2[Q^{z)-l]-39[nJ{z)-lf]/Q^{z). 

The dependence of the concentration on redshift and 
virial mass can be parametrized as 



= CA 



Mc 



(7) 



where the dimensionless parameters ca , cb , cc and the scale 
mass Mc depend on the cosmological parameters. For our 
fiducial model and the rest of this section, we use results 
from recent DM-only simulations of relaxed haloes (ca ~ 
9.23, CB = -0.09, cc = -0.69, Mc = 2 x lO^^fe-^Mo; Duffy 
et al. 2008; DOSfl which used WMAP5 cosmological pa- 
rameters, but we will consider other published relations in 
Section O 

Finally, one calculates the central density po by impos- 



ing 



p{r)4TVr dr = /gasMvir = (fb - /star)Mvir 



(8) 



where /gas and /star are the virial gas and stellar mass frac- 
tion, respectively. For /star, we will consider two limiting 
cases: 



5.7 X 10" 



Mv 



(9) 



5 X 1013 Mq^ 

following the work of Giodini et al. (20090, or /star — 0, 
depending on whether we want to quantify the effect (or lack 
thereof) on the mass budget of stars in low-mass systems. 
Since we are neglecting gas self-gravity, the choice of po (or, 
equivalently, of /star) does not affect the results on the M-T 
relation in Section [3. II However, it does affect the results on 
the L-T relation in Section [3.21 For the remainder of this 
section, we will assume /star = (i.e. /gas = /b). 

All parameters but one have now been fixed. The cen- 
tral pressure or, alternatively, the central temperature can 
in principle have any value. In reality, there is another con- 
straint that helps us fix Tq. 

A general attribute of the polytropic solution of the HE 
equation is that it is physically meaningful only between 
zero and a maximum radius rmax (Lea 1975). Beyond rmax, 
the temperature becomes negative, making the solution un- 
physical in that regiorQ. The existence of a maximum radius 
can be explained in the same way one explains the finite 

^ We assume Equation ((JJl to be valid for any 1 < /3 < 3 and for 
any < Mvir < lO^'^Mo. 

^ We approximate their results for the full sample of COSMOS 
selected groups only, extrapolating to higher masses, where the 
stellar fraction becomes negligible anyway, and assuming that the 
stellar fraction within r^oo is equal to the stellar fraction within 
rvir. This last point is an approximation, since we know that the 
gas mass increases with radius, but it does so only slightly. 
^ For some values of F (e.g. F = 3/2), the density is positive also 
for r > rmax, but the temperature becomes negative regardless 
of the value of the polytropic index. 
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truncation radius in classical polytropes with F > 6/5: a 
finite rmax means that the gravitational potential is strong 
enough to bind the gas in a finite region; an infinite rmax 
means that the pressure of the gas is high enough to make 
the gas either unbound or bound only in a limitless region. 
However, the analogy between the profiles in this work and 
those of classical polytropes cannot be taken too far: in clas- 
sical polytropes, r^ax = +00 if F < 6/5, regardless of other 
gas quantities; on the other hand, the polytropic solution 
given by Equation ((S]) or Q, for a given polytropic index, 
can have a finite or infinite maximum radius, depending on 
other gas quantities, as it is explained below. The difference 
arises from the different nature of 0: in this paper, the grav- 
itational potential is fixed and external, due to DM; in the 
classical polytrope case, is entirely due to gas self-gravity. 

By studying Equations ^ and ^ further, we note that 
rmax, once all the other parameters (virial mass, redshift, 
cosmological parameters, etc.) have been fixed, depends only 
on the gas temperature. The hotter the gas is, the larger 
Tmax becomes. In particular, by imposing the weak boundary 
condition Too > 0, there exists a threshold minimum value 
of the central temperature. 



ATvCScPcrsfJ-rnp F — 1 
fcB(/3- 1) 



1 /imp 



(10) 



above which the solution is physically meaningful at all radii, 
i.e. Tmax = -l-oo. This central threshold temperature, which 
depends on the system's gravitational potential (i.e. on its 
mass) and on the steepness of the polytropic relation, is a 
characteristic quantity of the system: To_thr is the minimum 
permitted central temperature, resulting from the require- 
ment of HE and non-negative temperature at all radii. As 
expected, a deeper gravitational potential requires a higher 
threshold central temperature. There exists also a maximum 
permitted central temperature, constrained by the boundary 
condition at large radii, usually given by the temperature 
of the intracluster medium or of the intercluster medium, 
depending on the type of system (e.g. group vs. cluster; iso- 
lated vs. non-isolated) under consideration. There is there- 
fore, for any given F, a relation between the mass of a sys- 
tem and a range of permitted central temperatures of its gas. 
This M~T relation will be fully exploited in Section[3]to con- 
strain gas and DM parameters. If To < To-thr, ''max is finite 
and can even be much smaller than the virial radius, depend- 
ing on the other parameters. If To — To-thr, ^max = +00 and 
all gas quantities at infinity vanish (and we recover the so- 
lution given in BIO). If To > To_thr, f'max = +00 and all gas 
quantities at infinity are positive. Imposing the boundary 
condition T(-l-oo) — Too, we have To = To-thr + Too- 

In Figure [T] we show the effect of changing To on the 
temperature radial profile, derived from Equations ((3| and 
((4| and from the polytropic relation, for a few cases of /3. 
The solution with To — To-thr is the only case for which the 
temperature becomes zero at infinity, and it does so with a 
logarithmic slope at large radii d(logT)/c?(logr) ~ — 1 for 
P > 2 and ~ — /3 + 1 for /3 < 2. This strong dependence 
on P will be useful in Section fS.l.Sl where we constrain the 
DM logarithmic slope at large radii ( = —fi — 1. Notice also 
that, in this particular case, the shape of the profile is in- 
dependent of the polytropic index. When To 7^ To_thr, the 
temperature profile T{r)/To depends on the polytropic in- 




1 2 3 

logio(r/rj 

Figure 1. Dependence of the polytropic solution to the HE equa- 
tion (Equations [3] and |4jl on the gas central temperature, for 
three different values of /3: 3/2 (upper, red bundle of curves), 2 
(NEW case; middle, black bundle of curves) and 5/2 [lower, blue 
bundle of curves). The black, vertical, dot-dashed line denotes 
the DM scale radius rg, whereas the vertical band denotes the 
virial radius, for a range of systems at redshift ^ = and with 
lO^^M© < Mvir < lO^'^M©. The temperature profile depends 
strongly on the value of /3. When To = To_thr [solid curves), the 
logarithmic slope at large radii d(logT)/d(logr) is ~ — 1 for P > 2 
and ~ — /3 -|- 1 for /? < 2. Additionally, at very large radii there 
is a strong dependence on To: the dotted [dashed) curves show 
the temperature profile when the central temperature is increased 
(decreased) by 1 per cent from To_thr- 



dex. However, when F varies, so do To_thr and its fractions 
(e.g. 0.99 To-thr and 1.01 To-thr, shown in the figure), there- 
fore Figure [1] looks exactly the same for any value of F > 1. 
If we slightly increase or decrease the central temperature, 
the effect on T(r)/To is quite dramatic at very large radii, 
whereas it is almost not observed within the virial radius. 

We can calculate the central temperature when one does 
not impose rmax = -l-oo, but imposes instead a finite rmax = 
Krvir, where k is some large integer. In this case, the new 
central temperature To-k is given by 



(1 + KCvi, 



\0~2 



To- 



1 - 



/3 — 2 KCvir(l + nCv 
ln(l -I- KCvir) 



'1/3-2 



ii 13 = 2. 



(11) 



It can be shown that the ratio To-^/To-thr > 0.9 for 
K > 10 (i.e. rmax > lOrvir), in the typical range 4 < Cvir < 10 
at 2: = 0, for a range of virial masses lO^'^ M© < Mvir < 10^^ 
Mq, for any 2 < /3 < 3. 

In extended systems, gas quantities at large radii are 
never really zero but are equal to those of the intracluster 
medium or of the intercluster medium, depending on the 
type of system one considers. Assuming global HE, we can 
compute the central temperature when one does not impose 
Too = 0, but imposes instead Too = Tq/F, where F is some 
large number. In this case, the new central temperature To_f 
is given by 



To- 



F 



To- 



F 



(12) 
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Figure 2. Dependence on rcdshift of several parameters, for a 
cluster of virial mass Mvir = 10^'' Mq, composed of NFW DM 
and gas with polytropic index F = 6/5 and central temperature 
^bTo = 2.5 keV. Top panel: we plot the temperature of the 
cluster To (dotted line), its central threshold temperature rQ_tijr 
(solid line) and its temperature at infinity Too (dot-dashed line). 
Notice how feB^O— thr does not depend very strongly on rcdshift, 
varying from ~ 2 keV to ~ 5 keV in the 2: = 0-3 range. Middle 
panel: wo plot the gas parameter Agas (solid line) in relation to 
the ratio r/(r— 1) (dotted line). Bottom panel: we plot the virial 
radius rvir (dotted line) and the maximum radius rmax (solid line) 
of the cluster. All panels: the dashed, vertical line denotes the 
threshold redshift ^thr S't which To = To_thri = 0, Agas = 
r/(r — 1) and the function rmax (2) has its vertical asymptote. 
When we impose Too = 0, 2: = Zthr is the only redshift at which 
the cluster can exist in HE at all radii. If we instead allow for 
any Too > 0, then the cluster can exist in HE at all radii in the 
2 = 0-2thr range. 



The ratio To_p/To_thr < 1-01 for F > 100 (i.e. Too < 
To/100), for any 1 < /3 < 3. In clusters of galaxies, where 
F >> 1, we can safely assume Too = 0. Groups of galaxies, 
however, can have a much smaller F, and the assumption 
of a zero temperature at infinity cannot be used. Written 
in a different way: in clusters of galaxies. Too is negligible 
with respect to To_thr (i.e. To = To-thr + Too ~ To-thr). In 
groups of galaxies, on the other hand, a non-zero Too cannot 
be neglected, and we will see in Section [3.1.21 its effects on 
the M-T relation. 

Considering quantities at infinity, like Too , is mathemat- 
ically convenient but obviously not realistic. For example, 
the HE assumption may break down at ^ r^h, where there 
may also be a virial shock. In that case, a more complete ap- 
proach might include the solution of the Rankine-Hugoniot 
jump conditions at the shock front, in order to obtain proper 
boundary conditions for the polytropic solution. A similar 
approach is required if one assumes the DM profile to have 
a sharp truncation radius at rvir (e.g. ©striker et al. 2005). 
However, note that while numerical simulations can show 
virial shocks at ~ rvir along particular directions from the 
cluster centre (e.g. pointing to filaments), when the cluster 
profiles are spherically averaged the simulations do not show 
strong DM or gas pressure and density discontinuities at rvir 
(e.g. Molnar et al. 2009). Moreover, given the uniqueness of 
the solution of the HE differential equation, once the Dirich- 



Figure 3. Dependence of the parameters of the cluster described 
in Figure |2] on its gas central temperature. The thick lines are 
the same lines shown in Figure [2] The arrows and the bands 
show what changes when we vary the central temperature from 
O.9T0 to 1.1 Tq. When increasing the central temperature, Agas 
decreases. Too and rmax increase and the threshold redshift 2tiir 
increases. Notice the existence of a minimum central temperature 
To— mini below which the cluster cannot exist in HE at all radii, 
at any redshift. 
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Figure 4. Dependence of the parameters of the cluster described 
in Figure [2] on its gas polytropic index. The thick lines are the 
same lines shown in Figure [2] The arrows and the bands show 
what changes when we vary the polytropic index from 0.98 F to 
1.02 r. When increasing the polytropic index, To_thr increases, 
Too and rmax decrease and the threshold rcdshift Ztlir decreases. 
Notice the existence of a maximum polytropic index Fmax, above 
which the cluster cannot exist in HE at all radii, at any redshift. 



let boundary condition has been enforced, the imposition of 
T(-l-oo) = Too can be mapped to an equivalent boundary 
condition T^ at some finite radius ^ rvir, with ^ > 0: 



To 



1 - (r - l)Agas [1 - ln(l eCvir)/eCvir] /F 

1 - (r - i)Agas/r 



(13) 
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for (3 = 2 and 



r 



/3-2 (€Cvir)(l + Cc 



(14) 



i-(r-i)Agas/r 

for /3 7^ 2. These formulas will be useful in Section [3.1.2l It 
is also worth remembering that observational quantities like 
M500 and T500 (see Section \2A\ do not depend strongly on 
the exact solution at > rvir, because of the low gas density 
and emissivity at large radii. 



2.2 Analysing the polytropic solution 

In this section, we show the behaviour of r-max and To_thr 
in a couple of different cases. In Figure [21 we show the de- 
pendence of Tmax as a functiou of redshift for a cluster of 
virial mass Mvir ~ lO^^M©, with gas of polytropic index 
r = 6/5 and central temperature kBTo = 2.5 keV, /gas = /b, 
and with DM described by an NFW profile and concentra- 
tion given by DOS. We also plot the virial radius and the 
threshold central temperature to allow for comparison. No- 
tice how the threshold central temperature slowly increases 
by only a factor of ~2 (from AiBTb-thr ~2 to ~5 keV) in the 
z = 0-3 rang4f|: for a fixed virial mass, the strength of the 
gravitational potential only slowly increases with redshift. 
Finally, we also show the gas parameter Agas, to compare 
it to r/(r — 1), and the temperature at infinity Too. The 
dashed, vertical line in all panels denotes the asymptote of 
the function rinax{z) and also denotes the threshold redshift 
2thr at which Agas = r/(r - 1), To = To_thr and = 0: 
below Zthr, the solution is physically meaningful at all radii. 
Above zthr, the radius rmax ^ +00 and decreases as redshift 
increases. Assuming HE at all radii, this threshold redshift 
is the only allowed redshift for this cluster, if we strictly im- 
pose Too =0, or is just a maximum redshift, if we allow for 
any non-negative Too. 

In Figures [3] and m we show what changes when we vary 
the central temperature (by ±10 per cent) or the polytropic 
index (by ±2 per cent) of the gas, respectively. If we de- 
crease the central temperature, the threshold redshift Zthr 
decreases, because the gas parameter Agas increases, and 
there is a minimum central temperature (in this example 
fcBro-min — 1.8 keV) uudcr which a cluster with a given 
virial mass and polytropic index cannot exist in HE at all 
radii, at any redshift. Analogously, if we increase the poly- 
tropic index, the threshold redshift decreases, because the 
threshold central temperature increases, and there is a maxi- 
mum polytropic index (in this example Fmax — 13/10) above 
which a cluster with a given virial mass and central temper- 
ature cannot exist in HE at all radii, at any redshift (see also 
Arieli & Rephaeli 2003). The existence of a minimum central 
temperature and of a maximum polytropic index is due to 
the fact that the threshold central temperature at redshift 
2 = is finite and is an increasing function of F. Notice that 
Fmax increases, for a fixed virial mass, if we increase the cen- 
tral temperature. However, there is an additional constraint 
on the polytropic index, given by the condition of convec- 
tive stability. The solution given in Equations ((3| and Q is 



° Note that, at high redshift, some of the assumptions of our 
model, including HE and isolation, are less valid. 



convectively unstable when {p/P)dP/dr < 'ydp/dr, which 
translates to F > 7, where the adiabatic exponent 7 = 5/3, 
for a completely ionized, ideal, monatomic, non-degenerate 
gas. 

On the other hand, if we increase To, the threshold 
redshift increases. The same happens if we decrease the 
polytropic index of the gas, because the threshold central 
temperature decreases. This trend can be seen also in an- 
other way. When we decrease the polytropic index towards 
its isothermal limit Fum = 1, the relation between pres- 
sure and density can be written as P{r) = [p{r)/po\To. 
In this limiting case, Equation Q is no longer valid for 
13 — 2, and the solution to the HE equation, using once 
again the boundary condition p(0) = po, become^ p(r) = 
po exp[-Agaa(l - ln(l -I- r/rs)/{r/ra))] (Makino, Sasaki & 
Suto 1998). It is clear from this solution that the gas quan- 
tities are never negative at any radius, therefore there is no 
minimum central temperature requirement: To-thr = 0. 

In the example shown in Figure [51 2thr ~ 0.85 for 
F = 6/5. We will see in Section (3.1.11 that the polytropic 
index has typical values between 6/5 and 13/10. This means 
that the threshold redshift computed for F = 6/5 is indeed 
a maximum value for the example shown, for any reason- 
able value of the polytropic index. If we could observe with 
enough accuracy the central temperature of the cluster gas, 
we could then obtain some useful information. Using the 
numbers above, if we for example observed a cluster of virial 
mass Mvir = lO^'' M© with a gas with polytropic index 
F > 6/5 and central temperature kuTo = 2.5 keV to be at 
redshift 2 > 1, then we could safely assume that the gas 
is not in HE at all radii. Alternatively, if we for example 
observed a cluster with the same virial mass and central 
temperature at redshift 2 ~ 0.5 and knew by some means 
that the gas is in global HE, then we could safely assume 
that it has a F > 6/5. Unfortunately, accurate observations 
of To are very difficult if not in many cases impossible, due 
to resolution limits of current X-ray telescopes. However, the 
same reasoning translates to the concept of average temper- 
ature, defined in Section [2.41 This will be evident in Section 
13.11 where we show the dependence of the M-T relation 
on several parameters, including the polytropic index and 
redshift. 

We now study the dependence of the HE-polytropic so- 
lution on the mass of the system. In Figure [5l we show 
the same parameters of Figure [2l this time as a function 
of virial mass, for a cluster at redshift 2 = with NFW 
DM, /gas ~ fh, and gas with polytropic index F = 6/5 and 
central temperature ksTo = 4 keV. Notice how the thresh- 
old central temperature increases significantly as the virial 
mass increases: for a fixed redshift, the strength of the grav- 
itational potential quickly increases with virial mass. As a 
consequence, rmax strongly depends on virial mass, changing 
from r-max ~ rVir to +CO in a very small mass range. 

In Figures[6land[7l we show what changes when we vary 
the central temperature (by ±10 per cent) or the polytropic 
index (by ±2 per cent) of the gas, respectively. The effect is 



^ In the /3 7^ 2 case, the isothermal solution is 



p{r) = po exp 



/3-2 (r/rs)(l + r/rs)'3-2 



> 0. 
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Figure 5. Dependence on virial mass of several parameters, for 
a cluster at redshift 2: = 0, composed of NFW DM and gas with 
polytropic index F = 6/5 and central temperature k^To = 4 keV. 
Top panel: we plot the temperature of the cluster To [dotted 
line), its central threshold temperature To_thr {solid line) and its 
temperature at infinity Too {dot-dashed line). Notice how To_thr 
depends strongly on virial mass, changing by more than one or- 
der of magnitude in the Mvir = 10^^ - 10^" M© range. Middle 
panel: we plot the gas parameter Agas {solid line) in relation 
to the ratio r/(r — 1) {dotted line). Bottom panel: we plot 
the virial radius r^i^ {dotted line) and the maximum radius rmax 
{solid line) of the cluster. All panels: the dashed, vertical line 
denotes the threshold virial mass Mvir— thr which Tq = To_thri 
Too = 0, AgaB = r/(r - 1) and the function '"maxC-Mvir) has its 
vertical asymptote. When we impose Too = 0, Mvir = Mvir— thr 
is the only mass that the cluster can have to exist in HE at all 
radii. If we instead allow for any Too > 0, then the cluster can 
exist in HE at all radii with a virial mass A/vir = 0-Mvir-thr- 

similar to what shown in Figures |3] and |4l except for the fact 
that there is no To-min or Fmax: for every small central tem- 
perature or large polytropic index, there is always a small 
enough virial mass for which the DM gravitational poten- 
tial is not strong enough to bind the gas in a finite region. 
Remember, however, that F < 7, because of the convective 
stability requirement. 

2.3 Comparison to previous studies 

In Section 12.11 we describe the structure of gas in clusters 
and groups of galaxies of known virial mass and redshift, 
by assuming HE at all radii between an ideal polytropic gas 
and DM. If the DM potential is known, the HE-polytropic 
solution is completely described once we fix its three free 
parameters: F, po and Tq. There are several ways to "close 
the problem", which are all mathematically equivalent. Here 
we briefly summarise the assumptions we made in Section 
12.11 and then contrast them with those made in prior work. 
This will help put into better context our results in Section 
[3] and those of other authors. 

For a given polytropic index, we obtain a relation be- 
tween the gravitational potential of a system and a range of 
permitted central temperatures of its gas, by simply requir- 
ing (i) spherical symmetry, (ii) HE between DM (described 
by a generalised NFW profile) and an ideal polytropic gas 



Figure 6. Dependence of the parameters of the cluster described 
in Figure [5] on its gas central temperature. The thick lines are 
the same lines shown in Figure \5\ The arrows and the bands 
show what changes when wc vary the central temperature from 
O.9T0 to 1.1 Tq. When increasing the central temperature, Agas 
decreases. Too and rmax increase and the threshold virial mass 
Afvir-thr increases. 
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Figure 7. Dependence of the parameters of the cluster described 
in Figure [5] on its gas polytropic index. The thick lines are the 
same lines shown in Figure [S] The arrows and the bands show 
what changes when we vary the polytropic index from 0.98 F to 
1.02 F. When increasing the polytropic index, To-thr increases, 
Too and rmax decrease and the threshold virial mass Mvir— thr 
decreases. 

at all radii, and by imposing (iii) that the temperature at 
large radii is non- negative (see Equation [TO] and related dis- 
cussion). 

KSOl (see also KS02) presented a similar description 
for the intracluster medium in clusters, by also assuming (i) 
and (ii), and additionally requiring (iv) that the gas density 
profile is equal to that of DM at large radii, (v) the DM large- 
radii logarithmic slope C, = —3 and (vi) the concentration is 
known, for a given mass and redshift. 

A06 (see also Ascasibar et al. 2003) also studied cluster 
models of polytropic gas, by once again requiring (i), (ii), (v) 
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and (vi), additionally (vii) taking the polytropic index as a 
given (from numerical simulations and observations), and 
imposing (viii) that the temperature at infinity is always 
zero, regardless of the system's mass or environment. 

Ostriker et al. (2005) (but see also Bode et al. 2007, 
2009) studied similar scenarios, by also assuming (i), (ii), 
(v), (vi) and (vii). Moreover, (ix) they constrained the cen- 
tral temperature by matching the external surface pressure 
to the momentum flux from the infalling gas at the virial 
radius. 



2.4 Average temperature and the issue of excision 



In this section, we connect the results given in Section[27T]to 
observations, by means of defining the average temperature. 
As shown in the previous sections, for a given polytropic 
index, we are able to obtain a relation between the gravi- 
tational potential of a galaxy group or cluster and its gas 
central temperature, once a boundary condition (e.g. Toe) 
has been fixed. Naturally, the central temperature is not an 
easy observable, because of resolution limits, but it is possi- 
ble to apply the same reasoning and results to the concept 
of average temperature. We computed the average temper- 
ature following the notation of KSOl. In general, one can 
write 



//^' w{r)T{r)4:Tvr^dr 
J^^ w{r)4-Kr'^dr ' 



(15) 



where the excision radius ro is either zero or a fraction of an 
overdensity radius (e.g. r^ir, ^200, r^oo, or other), the trunca- 
tion radius ri is some other fraction of, or equal to, an over- 
density radius, and li'(r) is the weight function. If the weight 
function w = 1, Tav = Tvwa is the volume- weighted average 
temperature; ii w — p, Tav ~ Tmwa is the mass-weighted 
average temperature; if to = mricAN, where rio is the num- 
ber density of electrons per unit volume, ni is the number 
density of ions per unit volume and An is the cooling func- 
tion, Tav = Tewa is the emission-weighted average tempera- 
ture. Although some authors use the mass- weighted average 
temperature, most observers use the emission-weighted av- 
erage temperature instead, usually approximating the cool- 
ing function as An oc p^T^^^. Accordingly, we also use the 
emission-weighted temperature here. However, since we also 
investigate the group regime, where temperatures are lower, 
we use a more general version of the cooling function, de- 
scribed in detail in Appendix|B] When using both definitions 
of the cooling function, we find the results to be basically 
identical at high temperature (as expected, since cooling is 
dominated by Bremsstrahlung emission) and only slightly 
different at low temperatures. 

The choice of ro and ri is an important one. The trun- 
cation radius ri is usually determined by observational con- 
straints. At large radii, both gas temperature and density are 
very low, therefore very difficult to quantify accurately. How- 
ever, exactly because temperature and density are very low, 
Tewa does not change appreciably for radii greater than rsoo , 
except for the highest mass systems, and does not change at 
all for radii greater than r^ii. Even when future telescopes 
(e.g. eROSITA; Predehl et al. 2007, 2010) will permit us 
to study the outer regions with better accuracy, a larger ri 
should not significantly change the average temperature. 



The choice of the excision radius ro is less straightfor- 
ward. Most studies justify excising the central region be- 
cause it is significantly affected by radiative cooling and by 
other physics of which we do not know the details, e.g. feed- 
back. The positive result of excision is the large decrease 
in scatter in the X-ray scaling relations, at all masses (e.g. 
Markevitch 1998), which is highly useful for cosmology mea- 
surements. However, excision also throws out potentially im- 
portant physical information relevant to understanding in 
detail the formation of the cluster, such as the effects of 
cooling, feedback, concentration, etc. Therefore, it is useful 
to study how the choice of ro affects the results. In Section 
13.1.41 we will see an example of how significant excision can 
be. 

Cognizant of all these considerations, in this paper 
we choose to use the integration limits of V09 (see also 
Kravtsov, Vikhlinin & Nagai 2006) and define the average 
temperature af°l 



T500 = 



JQ. 15 rsoo 



An T (r) 47rr 2 dr 



J 0.15 rsoo 



(16) 



When studying the effect of excision, we compute the 
average temperatures Tsoo-no-cxc and Tvir-no-exc, simply 
by changing the lower integration limit to ro — 0. 

The average temperature described here will be used in 
Secti0n r3.ll where we calculate the M-T relation for groups 
and clusters of galaxies and compare it to recent X-ray ob- 
servations. 



SCALING RELATIONS AND COMPARISON 
TO OBSERVATIONS 



3.1 The M-T relation 



In Section 12.11 we effectively obtained a relation between 
the temperature profile and the mass of clusters and groups 
of galaxies, once the temperature at large radii was fixed. 
For a proper comparison to observations, we defined the 
average temperature in Section [2.41 Using this quantity, we 
can then construct an M-T relation that may be compared 
to observations. 

In our model, the resulting M-T relation is a function 
of five parameters: the polytropic index F, the temperature 
at large radii T5 (or, equivalently. Too), the DM logarithmic 
slope at large radii ( (or, equivalently, /3), the concentra- 
tion Cvir and the redshift z. To develop intuition for how 
these parameters affect the M-T relation, in the following 
sub-sections we will fix four of these quantities and study 
the dependence of the M~T relation on the remaining fifth 
parameter. Comparing the predicted M-T relation to that 
from X-ray observations, which have improved significantly 
in recent years, leads to constraints on the allowable ranges 
of parameter values. A thorough examination of the allowed 
range would require a simultaneous variation of all parame- 
ters at once, to allow for possible degeneracies. Such a study 

When computing average temperatures within other overden- 
sity radii, wc keep the excision radius the same. For example, the 
integration limits for the computation of Tvir are ro = 0.15 ^500 
and ri = rvir- 
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is beyond the scope of this paper, although where we can, 
we will note the more obvious degeneracies. 

Before we embark on our study, some considerations are 
in order. We determine the temperature of a galaxy group 
or cluster by solving the thermally-supported HE equation. 
However, it is clear from simulations (e.g. Nagai, Vikhlinin 
& Kravtsov 2007) that random gas motions and rotation 
can contribute up to 10 per cent of the pressure support, 
leading to an over-estimation of the gas temperature needed 
for HE. On the other hand, the mass computed from X- 
ray observations is the thermally-supported HE mass (e.g. 
Fabricant, Lecar & Gorenstein 1980), resulting in an under- 
estimation of the true mass within rsoo of up to 10 per cent, 
for relaxed systems (e.g. Lau, Kravtsov & Nagai 2009; see 
also Mahdavi et al. 2008). Since the M-T relation is a power 
law, these two effects cancel each other and the comparison 
between our model and observations is meaningful (see also 
KSOl). The same reasoning applies to any other form of 
non-thermal support, including cosmic rays and/or magnetic 
fields (see e.g. Parrish et al. 2011). 

More importantly, in our model we assume a complete 
polytropic relation (i.e. the polytropic index is the same at 
all radii). Recent numerical studies (e.g. Shaw et al. 2010; 
Battaglia et al. 2011a,b) show instead a radial dependence. 
We note, however, that the variation of F with radius is not 
very large, especially in the range 0.15 rsoo -'"soo, where we 
calculate the average temperature (see Section [2.4|) . These 
same studies (Battaglia et al. 2011b) also show a slight de- 
pendence of F on redshift, and we will address this possibility 
in Section [3. 1.5 1 
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Figure 8. Dependence of the theoretical M-T relation on the 
polytropic index, for systems at redshift 2 = {E{z) = 1), with 
polytropic gas in HE with NEW DM (C = -3) at all radii, con- 
centration given by DOS and Too = 0. The relation can always 
be well approximated by a single power law. We plot the M—T 
relation for seven evenly-spaced values of 11/10 < F < 17/10. 
The black, solid lines are the results from our theoretical model, 
whereas the red, dotted lines show the empirical approximation 
from Equation 11171 . Notice how the polytropic index has a sig- 
nificant effect only on the normalisation. For a fixed mass, the 
increase of F results in a higher average temperature. In this and 
in the following figures, we assume our model scales self-similarly 
with redshift. We will see in Section |3. 1 . 51 that this assumption is 
well justified. 



3.1.1 Dependence on the polytropic index 

In this section, we consider systems at redshift z — 
{E{z) = 1), described by an NFW DM profile (C = -3), 
with concentration given by DOS and a polytropic gas with 
Too = 0, and study the dependence of the M-T relation on 
the polytropic index F. 

Figure [S] shows this dependence, when assuming that 
systems at all temperatures share the same value of F. We 
first notice that the relation can be very well approximated 
by a single power law, over more than three orders of mag- 
nitude in both mass and temperature (not shown in the 
figure). Moreover, the slope of such function does not no- 
tably change with F: the polytropic index has a significant 
effect only on the normalisation. With all these considera- 
tions, we can then approximate the local M-T relation with 
the fittin j"! formula 

i?(z)Msoo = A/„ (jT^)' (^)^ (17) 

where e = 6/5, a = 1.46, Mo = 3.9545 x lO" H'^Mq, 
ksTo — 5 keV and fceTsoo is in keV. This approximation 
closely reproduces the results of the model for typical values 

For this and some other fits in the paper, we performed 
a non-linear least-squares fitting with the IDL routine MPFIT 
(Markwardt 2009). Curiously enough, if we plot Mvir versus 
^500 > we obtain a relation with a self-similar slope: Mvir = 
Mo[r/(r - l)]^''5(T5aa/To)3/2, where Mq = 7.09 x 10^^ H-'^Mq 
and fcsTo = 5 keV. 



of 11/10 < F < 17/10. In this and in the following sections, 
we assume our model scales self-similarly with redshift. We 
will see in Section r3.1.5l that this assumption is well justified. 

The effect of varying the polytropic index is consistent 
with what was already shown in Figures |4] and [7] for a given 
virial mass, when F increases, the threshold central tempera- 
ture increases, and so does the average temperature. When F 
decreases, the average temperature quickly approaches zero 
proportionally to ~[(F — 1)/T]'*^^, as the polytropic relation 
tends to its isothermal limit (i.e. F = 1). Notice that, be- 
cause of the dependence on the F/(F — 1) term, it is much 
easier to discriminate between lower values of the polytropic 
index than between higher values. 

We now compare the theoretical M-T relation to re- 
cent X-ray observations of galaxy clusters and groups. The 
red, dot-dashed line in Figure [5] shows the best-fitting M- 
T relation calculated by V09 for a sample of seventeen, 
low-redshift, relaxed clusters: E{z)M5oo = (3.02 ± 0.11) x 
lO"/i-iM0(fcBr5oo/5 keV)i '^3±°-°*, valid over the temper- 
ature range ksT ~ 1.5-12 keV. Since we do not have ac- 
cess to the exact measurement and error analysis of V09, 
we cannot make a definitive statement on the range of F 
consistent with V09. Assuming zero covariance between the 
magnitude and the slope of the published M-T relation, we 
find that a model with a mass-independent polytropic index 
1.22 < F < 1.24 is consistent with the data. Without this 
assumption, we conservatively estimate 6/5 < F < 13/10. 
Such values of F are also consistent with the results ob- 
tained by both observations and simulations (e.g. Bode et 
al. 2009 and references therein). 
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Figure 9. Constraints on the polytropic index, for systems at red- 
shift z = {E{z) = 1), with polytropic gas in HE with NFW DM 
(C^ = —3) at all radii, concentration given by DOS and Too = 0, 
via comparison between our theoretical model and observations 
by V09 and S09 {red, dotted and turquoise, long-dashed line, re- 
spectively). We plot the theoretical M-T relation for F = 6/5 
{upper, black, solid line) and F = 13/10 {lower, blue, dot-dashed 
line). A model with a mass-independent F with a value between 
1.22 and 1.24 is consistent with the l-cr errors of V09. No model 
with a constant F is consistent with S09. For consistency with 
S09, with all other parameters fixed, we need to assume a mass- 
dependent polytropic index (see Eguation llSp . 



Observations of low-mass systems show increasing scat- 
ter and a steepening of the M-T relation. It is not clear 
if the M-T relation for all (low- and high-mass) systems 
should be described by a single power law (only steeper 
than the cluster-only case), or if there is a real break in the 
relation at low masses. The existence of such break is still 
under debate (e.g. Eckmiller, Hudson & Reiprich 2011). Sev- 
eral authors (e.g. Nevalainen, Markevitch & Forman 2000; 
Finoguenov, Reiprich & Bohringer 2001; Sanderson et al. 
2003) suggest its presence at around k^T ~ 1 keV, but 
other authors (e.g. S09) do not find conclusive evidence for 
it. We will consider the possibility of a broken power law 
in Section 13.1.21 and focus here on comparing our model 
to a single, steeper power law. The turquoise, long-dashed 
line in Figure |9] shows the best-fitting M-T relation calcu- 
lated by S09 for a sample of thirty-seven, low-redshift, re- 
laxed groups and clusters (including fourteen clusters from 
V09): E{z)Mtmo = (1.26 ± 0.07) x 10"/i-^MQ(fcBr5oo/3 
keV)^'^^^"'"*, valid over the temperature range kBT ~ 0.7- 
12 keV. In this case, our model with a constant polytropic 
index is clearly not consistent with the data. 

A possible explanation for this discrepancy is that the 
polytropic index is not the same for all systems, but instead 
depends on the mass of the system. By direct comparison 
of the theoretical model to the best-fitting M-T relation 
by S09, we obtain an empirical relation between F and the 
virial mass of local groups and clusters of galaxies, 

F(MviO = 1.305 - 0.073 logio (j^^) - (18) 



which is valid over the combined V09 and S09 temperature 
ranges - k^T ^ 0.7-12 keV - and where Mvir is in M0. 

The polytropic index slowly decreases with mass. For a 
proper comparison with previous literature, we convert this 
equation to a relation between F and concentration, by using 
Equation 0: 

r(cviO = 1.305 + 0.816 logio (H) , (19) 

where 6.8 is the concentration of a local cluster of virial mass 
equal to 10^^ M©, according to D08. 

KSOl (see also KS02) and A06 each obtain a relation 
between F and Cvir- Although in all three relations the poly- 
tropic index increases with concentration, the slope and nor- 
malisation are quite different. A06, by using data from their 
own simulations, find a polytropic index which is almost con- 
stant with concentration: F = 1.145 -I- 0.0036 Cvir. KS02, by 
requiring that the gas profile at large radii is equal to that 
of DM, find a slightly steeper slope than A06, obtaining the 
relation F = 1.137 -f 0.0894 ln(cvir/5) - 0.00368(cvir - 5), 
valic0 for 1 < Cvir < 25. In both cases, the polytropic 
index in the combined temperature range of S09 and V09 
(/cbT ~ 0.7-12 keV, Cvir ~ 5-7.5) is almost constant and 
close to 1.15, a value slightly lower than that of current ob- 
servations and numerical simulations (e.g. Bode et al. 2009 
and references therein). In that same temperature range, our 
relation is steeper than both A06 and KS02 (and KSOl), 
with the polytropic index varying from ~13/10 in the group 
regime to ~6/5 in the cluster regime. If we assume our the- 
oretical model and the observations we use as comparison 
(S09) to be correct, one implication of this discrepancy with 
KSOl and KS02 is the following; although qualitatively rea- 
sonable, the assumption that the slopes of gas and DM pro- 
files at large radii are equal cannot be used to construct an 
accurate model of gas in groups and clusters of galaxies. 

It is not immediately clear why the polytropic index 
should depend on the mass (or, equivalently, concentra- 
tion) of the system. Following Bertschinger (1985), F is 
determined by events like cluster collapse and relaxation, 
which might depend on mass scale and redshift, so that 
such a result may not be surprising. Also, feedback efficiency 
might be a strong function of halo mass and cause a mass- 
dependent polytropic index. However, we remind the reader 
that this constraint on F was obtained after fixing all other 
free parameters. In Section 13.1.21 we present a different ex- 
planation for the steepening/breaking of the M-T relation, 
wherein a mass-dependent polytropic index is not necessary 
and the shape of the power law changes with the choice of 
the temperature at large radii. Of course, both these causes 
could be at work at the same time. 



3.1.2 Dependence on the gas temperature at large radii 

In this section, we consider systems at redshift z — 
{E{z) = 1), described by an NFW DM profile (C = -3), 
with concentration given by D08 and a polytropic gas with 



KSOl provide a simpler relation, F = 1.15 -I- 0.01(cvir — 6.5), 
valid only for the temperature range we are considering in this 
paper. 
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Figure 10. Constraints on the gas temperature at large radii, 
for systems at redshift 2 = {E{z) = 1), with polytropic gas in 
HE with NFW DM (C = —3) at all radii, concentration given by 
DOS and F = 6/5, via comparison between our theoretical model 
and observations by V09 {red, dotted line) and by S09 {turquoise, 
long-dashed line and data points with l-cr errors). We plot the 
theoretical M—T relation for Too = {black, solid line), Too = 
lO^K {red, short-dashed line), T^o = 5 X lO^K {green, dot-dashed 
line) and Too = lO^K {blue, triple- dot- dashed line). A model with 
Tex, = is consistent with V09 but not with S09. For consistency 
with S09, with all other parameters fixed, we need to assume 10® 
K < Too < 10^ K, in agreement with numerical simulations and 
observations of the warm-hot intergalactic medium. 

r = 6/5, and study the dependence of the M-T relation on 
the temperature at large radii. 

In Section [2] after imposing HE at all radii, we found 
a range of permitted values for the central gas tempera- 
ture, To — min To — max 5 where To — mill = To_thr is constrained 
by requiring Too = and To-max is given by some different 
boundary condition at large radii (Too 7^ 0). Here, we explore 
this range of permitted temperatures further, by studying 
the dependence of the M-T relation on T^, the gas temper- 
ature at a given radius ^rvir. As demonstrated in Section 
12.11 there is a one-to-one relation between T5 and Too, for 
any given value of ^. Therefore, we can perform all calcu- 
lations using Too as the boundary condition for simplicity, 
and then translate the results back to T5. 

As shown in Figure [TO] when imposing Too 7^ 0, the M- 
T relation ceases to be a single power law. Although there 
is basically no difference at high masses (for any reasonable 
value of Too), the M~T relation steepens at low masses, re- 
producing an effective break. The larger Too is, the higher is 
the mass (and temperature) at which the turn-off should be 
observed. A steepening/breaking in the M-T relation has al- 
ready been mentioned in Section TS.l.ll as a possible feature 
of the scaling relation of low-mass systems (e.g. Eckmiller et 
al. 2011). In Section [3.1.11 we considered the single-power- 
law description and proposed that a mass-dependent poly- 
tropic index could explain the M-T relation that fits both 
low- and high-mass quantities (e.g. S09). Here, we consider 
the broken-power-law description and explain the existence 
of the break with a non-zero Too. We note, however, that 
the explanations we give are independent of the descrip- 
tion of the M-T relation; a mass-dependent polytropic in- 



dex can reproduce a broken power law, whereas a non-zero 
Too can also explain a single, steep power law. Likely, these 
two causes are concurrent. 

We propose here two interpretations for the tempera- 
ture at large radii. If we assume Too to be the same for every 
system - as one would expect, if one thinks of all isolated 
groups and clusters to be embedded in the same interclus- 
ter medium - we can constrain the lower and upper limits 
of the temperature of the external medium, since we actu- 
ally observe a possible turn-off at around fceT ~ 1 keV (e.g. 
Sanderson et al. 2003; see also Dos Santos & Dore 2002). 
In order to match the observations of galaxy groups. Too 
must be between 10^ K and 10^ K. This in turn results in 
different values of Tj for different values of virial mass (and 
of ^), from using Equation (fT5)) (or Equation [Til if /? / 2). 
For example, if we set ^ = 5 (i.e. if we are interested in 
the temperature at five virial radii) and Too = 5 x 10® K, 
we obtain Tj ~ 0.5, ~ 0.8 and ~ 2 keV, for local systems 
with M500 = 10^^, 10" and lO^'^ /I'^Msoo, respectively. We 
note that the values for the temperature of the intercluster 
medium obtained in this way are in agreement with numer- 
ical simulations (lO'^-lO'' K; see e.g. Cen & Ostriker 1999; 
Dave et al. 2001) and observations (e.g. Williams et al. 2007) 
of the warm-hot intergalactic medium. 

Alternatively, we can relax the assumption of isolation 
and consider the environment of groups and clusters, effec- 
tively considering different values of Too for each system. 
Whereas local clusters, given their size, can almost always 
be considered as isolated systems, most galaxy groups, on 
the other hand, are believed to exist embedded in proto- 
clusters (or even in clusters). Therefore, the effective Too 
of groups tends to be higher than that of clusters, caus- 
ing a steepening of the M-T relation. Fossil groups, on the 
other hand, can be considered to be fairly isolated, hence 
to have a lower effective Too than non-isolated groups (cf. 
Khosroshahi, Ponman & Jones 2007). This range in the en- 
vironment of low-mass systems - from isolated, to embedded 
in proto-clusters, to embedded in clusters - could be one of 
the possible causes for the larger scatter of the M-T relation 
at low temperatures than at high temperatures (see Sec- 
tion [31113] for another possible cause). Likely, these two in- 
terpretations (mass-dependent and environment-dependent) 
co-exist. 

Finally, if we allow for the truncation of the DM gravi- 
tational potential or for the breaking of the HE assumption 
at a given, finite radius (e.g. ^ Tvir), the polytropic solution 
is not global any more. In this case, T^ can be interpreted 
as being related to the external surface pressure term in the 
virial theorem, as Ostriker et al. (2005) assumed in the case 

C = i- 

3.1.3 Dependence on the DM logarithmic slope at large 
radii 

In this section, we consider systems at redshift z = 
{E{z) — 1), described by a generalised NFW DM profile 
(BIO), with concentration given by DOS and a polytropic 
gas with a given F and Too = 0, and study the dependence 
of the M-T relation on the DM logarithmic slope at large 
radii 1^. 

To achieve this, we study the effects of varying /3 = 
— 1^ — 1 in Equation ([T]). The BIO profile is itself a simpli- 
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Figure 11. Constraints on the DM logarithmic slope at large 
radii, for systems at redshift z = {E{z) = 1), with BIO DM, 
concentration given by DOS and a polytropic gas with Too = 0, 
via comparison between our theoretical model and observations 
by V09 and S09 {red, dotted and turquoise, long-dashed lines, 
respectively). We plot the theoretical M—T relation for two poly- 
tropic models — F = 6/5 [black lines) and T = 13/10 {blue lines) 
— and three DM models - C, = —5/2 [dotted lines), f = —3 
(NFW; solid lines) and = —7/2 [dashed lines). Assuming a 
mass-independent polytropic index 6/5 < F < 13/10 for high- 
mass systems, we obtain —3.5 < C $ —2.7: the DM in galaxy 
clusters either follows the NFW profile or is slightly steeper at 
large radii. 



fied version of the general profile given in Suto et al. (1998). 
However, Suto et al. (1998) themselves, and several authors 
thereafter (including KSOl and KS02), use a simplified ver- 
sion of their own profile (see also Zhao 1996), by fixing the 
large-radii logarithmic slope to ^ = —3. The advantage of 
using the BIO profile is that we can compare our model to 
observations and obtain constraints on how DM behaves at 
large radii. See Appendix|X]for the computation of different 
overdensity masses and radii in the case of a BIO profile. 

Comparing the results of our model to recent X-ray ob- 
servations of clusters (V09), we can discriminate between 
different values of i^, assumed to be the same for systems 
at all masses. In Figure 1111 we show the M-T relation for 
two values of the polytropic index and for three values of 
Assuming that a mass-independent polytropic index has 
a value 6/5 < F < 13/10, we obtain -3.5 < C ^ -2-7: the 
DM in galaxy clusters either follows the NFW profile or is 
slightly steeper at large radii. The comparison to low-mass 
systems (S09) is less straightforward: a constant value of ( 
is not consistent with the data. However, more causes could 
be at work, including a mass-dependent polytropic index 
and/or a non-zero temperature at large radii (see Sections 
I3.1.1l and l3.1.2p . hence complicating the determination of ^. 

As expected, for a fixed mass, the temperature is higher 
for lower values of /? (less negative values of Q . This is easily 
understood by investigating Equation (|10|l or the definition 
of (t>o, which is proportional to l/(/3 — 1): the lower /3 is, 
the deeper the central gravitational potential and therefore 
the higher the central temperature has to be, in order for 
the gas to be unbound or bound only in a limitless region. 
The relation between the central temperature and the av- 



erage temperature is not, however, a linear relation. The 
parameter /3 also affects the shape of the temperature pro- 
file (through the gas parameter Agas , given by Equation [SJ . 
For this reason, the M—T relation does not shift indefinitely 
towards lower values of temperature (for a given mass). In- 
stead, for values of /3 > 5/2, it slightly shifts towards higher 
values of temperature, especially at lower masses (not shown 
in the figure). 

One can speculate on the interdependence between C, 
and r, hinted at by Figure 1111 Does the gas behave differ- 
ently depending on the gravitational potential - i.e. do we 
have a T[C,) relation? - or does the DM respond somehow 
to the gas - i.e. C(r)? Even though it is possible to think 
of scenarios where DM is significantly affected by baryons 
(e.g. adiabatic contraction being one of them), it is unlikely 
that such physics is this effective at large scales. If there is 
indeed a direct relation between the polytropic index of the 
gas and the logarithmic slope of DM at large radii, it is more 
likely that the gas is affected by DM: r(("). This is consis- 
tent with the explanation given by Bertschinger (1985) for 
the existence of a polytropic relation in the first place: gas 
adjusts itself according to the formation process of the DM 
halo, and the formation process has an effect on the value 
of C (e.g. Lu et al. 2006). 

3.1.^ Dependence on the mass- concentration relation 

In this section, we consider systems at redshift 2 = 
{E[z) = 0), described by an NFW DM profile (C = -3), 
and a polytropic gas with F = 6/5 and Too = 0, and study 
the dependence of the M-T relation on the concentration 

Cvir . 

Numerical simulations and observations show that con- 
centration is a function of mass and redshift (see e.g. Equa- 
tion [7|). We focus here on the dependence of concentration 
on mass. Given the vast amount of published results and 
the discordance amongst them, we chose three representa- 
tive mass-concentration relations with three fairly different 
dependences on mass: cb ~ —0.09 (DOS), cb ~ —0.4 (Okabe 
et al. 2010: OlO) and cb = -0.59 (Oguri et al. 2011: OH). 
Amongst these relations, there is substantial agreement in 
the value of the concentration of high-mass systems. Concen- 
trations for low-mass systems, however, differ significantl){3- 

In Figure 1121 we show the M~T relation for the three 
mentioned mass-concentration relations. Surprisingly, very 
different concentrations do not translate into significant dif- 
ferences in the Afsoo-Tsoo relation, except in the very high- 
temperature regime, which is difficult to analyse, since very 
few clusters live there. Different concentrations imply var- 
ied gravitational potentials, and one would expect the gas 
to respond accordingly, giving rise to sizeably different M- 
T relations. This is a clear example where excision discards 
important information. If we compute the Msoo-Tsoo-no-exc 

^•^ As in the previous sections, we assume Equation ((7]| to be 
valid for all virial masses between lO^^M© and IQ^^Mq. DOS used 
redshift-dependent DM-only simulations, whereas OlO and Oil 
used lensing observations of clusters at mean redshift Zmoan = 
0.23 and 0.45, respectively. Assuming the same redshift depen- 
dence for all three relations (i.e. cq = —0.69), we then normalised 
the results of OlO and Oil by dividing their parameter c\ by 

/dm(-'- ~^ ^mcan)^^. 
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Figure 12. Dependence of the theoretical M-T relation on the 
Cvir— Mvir relation, for systems at redshift ^ = (E{z) = 1), with 
polytropic gas in HE with NFW DM (C = -3) at all radii, V = 6/5 
and Too = 0. Black, solid line: concentration given by DOS. Red, 
dashed line: concentration given by OlO. Blue, dotted line: con- 
centration given by Oil. Top panel: excision (r'o = O.lSrsoo in 
the computation of T500 ) • Different mass-concentration relations 
do not cause significant differences in the M-T relation. Bottom 
panel: no excision (ro = 0). The effect is much stronger: excision 
clearly discards information. 



relation instead (by replacing ro — 0.15 rsoo with ro = 
in Equation I16p . we obtain a much stronger effect: at low 
masses, the temperature is higher for lower (more negative) 
Cb . The effect of concentration is sizeable only in the inner- 
most region of the system. By performing excision to remove 
scatter, one also removes any information on the concentra- 
tion. Moreover, part of the scatter of the M-T relation at 
low temperatures might be due to the intrinsic scatter in the 
mass-concentration relation. We additionally note that exci- 
sion also removes information on different DM inner slopes 
(e.g. Moore et al. 1998) or on the effect of stars on the grav- 
itational potential, especially at low masses (e.g. Capelo et 
al. 2010). 

Incidentally, we notice that the mass-concentration rela- 
tion might be partially responsible for the steepening of the 
M-T relation (with excision) at low masses. However, as we 
will see in Section [3.21 mass-concentration relations with a 
strong mass dependence (like OlO and Oil) are ruled out 
by the comparison of the model to L-T observations. 
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Figure 13. Dependence of the theoretical M-T relation on red- 
shift, for systems with polytropic gas in HE with NFW DM 
(<^ = —3) at all radii, concentration given by DOS, F = 6/5 and 
Too = 0. We plot the M-T relation for redshifts z = {black, 
solid line), 0.5 (turquoise, dotted line) and 1 (blue, dotted line). 
The overlap of these three lines shows that the scaling of the 
model with redshift is very close to being self-similar. We also 
plot the M-T relation at redshift z = I for two different cases 
of non-self-similarity (see text for more details): = —0.5 [red, 
dashed line) and 0.5 (red, dot-dashed line). These two relations 
are well fit by our model when F varies with redshift according to 
Equation II26I I. when rj is given by Equation II30I I {blue lines). The 
same (not shown in the figure) can be achieved using Equations 
1281 and 1311 . 



3.1.5 Dependence on redshift 

In this section, we consider systems described by an NFW 
DM profile (C ~ —3), with concentration given by DOS and 
a polytropic gas with F = 6/5 and Too = 0, and study the 
dependence of the M-T relation on redshift. 

If we assume all galaxy group and cluster potentials 
to have the same shape (i.e. a fixed concentration) then, 
from the virial theorem, oc Mac/^'Ac oc M^^^p5rK?an, 

where pmoan = AcPc is the mean density within rA^- Since 
Pc oc E^(z), we expect 



mZ'e'^H^) 

^ Ml/'Al(!{z)E'^'{z) 



if Ac 
if Ac 



fixed number, 
Avir(z). 



(20) 



If we assume a flat, matter-dominated Universe (i.e. 
f2o = 1), then E{z) oc (1 -I- z)^^^ and we recover the results 
by Kaiser (1986), where Ta, oc M'^^{1 + z). 

In Figure 1131 we compare the results of our model at 
high redshift {z = 0.5 and 2 = 1) to those at 2; = 0, finding 
that the scaling of the model with redshift is very close to 
self-similar, according to Equation (|20p . even when using a 
concentration that varies with mass and redshift. Qualita- 
tively, the effect of varying redshift is consistent with what 
already shown in Figure O for a fixed F, when z increases, 
the threshold central temperature increases, and so does the 
average temperature. We do not extend the model to higher 
redshift, because the assumptions of HE and isolation be- 
come less and less valid with increasing z, due to the fact 
that clusters are actively assembling at these epochs. 
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In the past decade, several high-redshift X-ray obser- 
vations of clusters have been conducted (e.g. Reichert et 
al. 2011 and references therein). In some cases, the self- 
similarity of the M-T relation has been observed (e.g. Ettori 
et al. 2004; Reichert et al. 2011), and the theoretical model 
with the parameters listed above is in good agreement with 
it. 

However, there exist observations (e.g. Jee et al. 2011; 
Clerc et al. 2011; see also numerical simulations by Short 
et al. 2010) which show a discrepancy in the normalisation 
of the M-T relation between the observed values and the 
values expected from self-similarity. 

If we assume that our model is accurate also at high red- 
shift, there are two plausible ways to recover the observed 
non-self-similar relations. Assuming that the polytropic in- 
dex does not vary with time, one needs the DM profile to 
depend on redshift (e.g. ^{z)). Assuming instead that the 
DM profile does not evolve, one needs some variation of the 
polytropic index with redshift (e.g. T{z)). These two causes 
might be concurrent. Here, we focus on the latter possibil- 
ity and present a simple method to obtain the F-z relation 
from observations. 

It is customary to parametrize the observed discrepancy 
from self-similarity by adding a redshift-dependent term to 
the M-T relation: 



E{z)MA,ocT^jM-Tiz), 



(21) 



where the value of a, assumed to be the same at all redshifts, 
does not affect the results of this section, and fM-T{z) varies 
according to the author (with /m-t(« = 0) = 1). 

We now parametrize the time-dependence of the poly- 
tropic index accordingly and write 



r{z)^rofriz), 



(22) 



where To = r{z = 0) is the local value of the polytropic 
index and fr{z = 0) = 1. 

By combining together Equations (fTf)) . (pi) and (f22)) . 
our theoretical model reproduces the discrepancy from self- 
similarity when 



friz) = 



ro/M-T(z)i/7(ro-i) 



(23) 



To [ro/M-T(z)i/7(ro-i)]-i' 

The most common parametrizations for the non-self- 
similar M-T relation are (Case A; e.g. Short et al. 2010; 
Clerc et al. 2011) 



E(z)Ma, xT^^{l + zf 
and (Case B; e.g. Reichert et al. 2011) 

E{zfMA, ocTK,, 

where ip and v quantify the discrepancy from the self-similar 
cas For self-similarity, one needs tp = and v = I. 



(24) 



(25) 



^"^ These two parametrizations give very similar results only in 
the range ^ = - 1, which is exactly the range we consider in this 
section. At ^ = 1, Equation l|24| l with i/i = 0.5 or —0.5 is similar 
to Equation I I25I I with v = 0.33 or 1.67, respectively. 



The above descriptions imply /m-t(2) = (1 -I- z)'^ and 
/m-t(z) = E{zy~", respectively, in Equation (pT|) . We 
now parametrize the time-dependence of the polytropic in- 
dex accordingly. We assume fr{z) — (1 + z)^ when us- 
ing /m-t(2) = (1 -f z)^ , and friz) — E{zY when using 
/m-t(2) = E{zY^^ , where rj and 9 quantify the evolution 
of r. For self-similar models, one has r] = 6 = Q. 
The resulting V -z relation is then (Case A) 



r(z) = ro(i + 2)7 

where 

'7logio(l + z) = logio 
or (Case B) 

r(^) = ro£;(^)'', 

where 



(i + 2)'*/7(ro-i) 



[ro(i + 2)'^/7(ro - 1)] - 1 



'login £( 



logi 



E{z) 



/(ro - 1) 



[roi5(2)(i--)/7(ro - 1)] ■ 



(26) 



(27) 



(28) 



(29) 



These relations are very weakly dependent on redshift, 
especially in the redshift range z = - 1 and for low values 
of ^. Hence, after fixing the other parameters (e.g. e = Fq = 
6/5, see Section |3.1.1|) . we can provide a relation between 77 

and -0, 



r; = -0.17V + 0.06V , 
and between 9 and v, 

6 = -0.126-f- 0.082 u-f 0.044 



(30) 



(31) 



Observations of non-self-similarity are still giving very 
disparate results (e.g. Reichert et al. 2011; Clerc et al. 2011). 
In the future, when better observations will be at hand and 
the value of ip will be known to higher accuracy than now, 
it will be possible to infer the rate of evolution of the poly- 
tropic index, by using e.g. Equations H26|) and H30|) . This will 
also let us compare our results to numerical simulation (e.g. 
Battaglia et al. 2011b). 

3.2 The L-T relation 

In this section, we derive the L-T relation from our model 
and compare it to the observed relation of a sample of re- 
laxed clusters: S(z)-^L500 = (6.71 ±0.34) x lO"" g cm^ s"^ 
(fcBr5oo/6 keV)^-^^*" " (Mil). To be consistent with how 
we defined the temperature (in Section [2. 4|l . we set 



L500 = / ncniAN47rr dr, (32) 

where the cooling function An is described in Appendix [B] 
In Section [37ll we determined the temperature of groups 
and clusters of galaxies by imposing HE at all radii, between 
a polytropic gas and a given DM gravitational potential. 
For that calculation, only two constraints were necessary. 
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Figure 14. Dependence of the theoretical L-T relation on the 
polytropic index and on the gas fraction, for systems at redshift 
2 = (E{z) = 1), with polytropic gas in HE with NFW DM (C = 
—3) at all radii, concentration given by DOS and Too = 0. Black, 
solid lines: T = 6/5. Black, dotted lines: V = 13/10. Red, dashed 
lines: observations by Mil. Top panel: /gas = /b (/star = 0). 
Bottom panel: /gas varies (/star is given by Equation |9ll . 

namely on the polytropic index and on the central gas tem- 
perature, because we did not include gas self-gravity in the 
calculations. In this section, when deriving the luminosity, 
a third constraint (on the central gas density) is necessary. 
Different authors use different methods: KSOl (and KS02) 
impose the gas density at large radii to be /b times the 
density of DM; Ascasibar et al. (2003) require the baryon 
fraction to never exceed the cosmic value; A06 assume the 
baryon fraction at three DM scale radii to be equal to Q.\^/Q.o; 
Ostriker et al. (2005; see also Bode et al. 2007, 2009) impose 
conservation of energy. 

In this paper, the constraint on po is given by imposing 
the baryon fraction at the virial radius to be equal to the 
cosmic value, independent of mass or redshift. Therefore, the 
model depends also on the stellar content, and we compare 
results assuming the stellar mass fraction is given either by 
Equation Q or by /star = (i.e. /gas = /b). 

In Figure 1141 we show the theoretical L~T relation 
for two values of the polytropic index and one mass- 
concentration relation (DOS), and compare it to recent ob- 
servations by Mil, for both cases of constant and mass- 
dependent gas mass fraction. We note how the L-T relation 
is well approximated by a power law only at high temper- 
atures, when assuming a constant gas mass fraction, and is 



Figure 15. Dependence of the theoretical L-T relation on the 
polytropic index, the gas fraction and the concentration, for sys- 
tems at redshift 2 = {E{z) = 1), with polytropic gas in HE with 
NFW DM (C = -3) at all radii and T^o = 0. Solid lines: F = 6/5. 
Dotted lines: F = 13/10. Blue lines: concentration given by OlO. 
Turquoise lines: concentration given by Oil. Dashed, red lines: 
observations by Mil. Top panel: /gas = /b {/star = 0). Bottom 
panel: /gas varies (/star is given by Equation [9ll . 

well described by a power law at all temperatures, when as- 
suming a mass-dependent gas mass fraction. The effect of a 
gas mass fraction dependent on the total mass of the system 
is clearly visible, especially at low masses. This is a poten- 
tial source of scatter, since different systems having the same 
total (low) virial mass might have slightly different stellar 
masses. 

When assuming a mass-dependent gas mass fraction, 
the L-T relation with F = 6/5 matches X-ray observations 
better than the case with F = 13/10. 

Focusing on other mass-concentration relations (e.g. 
OlO or Oil, see Figure [T^ . we immediately notice that the 
model does not recover the observed L-T relation by Mil. 
This might imply that concentration cannot depend strongly 
on mass. 



4 SUMMARY AND CONCLUSIONS 

We present a simple and flexible model for the gas and DM 
content of groups and clusters of galaxies, by assuming an 
ideal polytropic gas in HE at all radii with a DM gravi- 
tational potential, described by a generalised NFW profile. 
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Exploring the physical properties of the polytropic solution 
of the HE equation, we report on methods to constrain the 
three gas parameters (F, To and po) and two DM parame- 
ters ((^ and Cvir) of such a solution. With no prior knowledge 
of the polytropic index or the boundary conditions, we de- 
rive the scaling relations for galaxy groups and clusters from 
these first principles and compare them to recent X-ray ob- 
servations. This allows us to obtain powerful constraints on 
gas and DM parameters. In particular, 

• for systems at redshift z = 0, with NFW DM, concen- 
tration given by DOS and Too ~ 0, we derive a theoretical 
M-T relation with a logarithmic slope of 1.46 and normali- 
sation that depends rather simply on the polytropic index: 
[r/(r - (see Equation [nil; 

• comparing with recent X-ray observations of clusters of 
galaxies, we find that a model with a constant polytropic 
index with a value between 6/5 and 13/10 is consistent with 
the data; 

• further, when including data from low-mass systems, 
the assumption of a constant polytropic index is not con- 
sistent with the observations. A possible explanation for the 
discrepancy is a mass dependence in F such that it decreases 
with mass (and increases with concentration), as derived in 
Equations HH]) and (HSJ; 

• for systems at redshift 2 = 0, with NFW DM, con- 
centration given by DOS and F — 6/5, we derive a the- 
oretical M-T relation that is a broken power law, where 
the break depends on the value of the temperature at large 
radii. Comparing this case to X-ray observations of groups, 
we constrain the temperature at large radii to be 10^ K 
^ Too ^ lO'' K, consistent with numerical simulations and 
observations of the warm-hot intergalactic medium; 

• for systems at redshift z — 0, with a generalised NFW 
DM, concentration given by DOS and Too = 0, we com- 
pare the theoretically derived M-T relation to recent X- 
ray observations of clusters of galaxies. Assuming a mass- 
independent polytropic index 6/5 < F < 13/10 for high- 
mass systems, we find that —3.5 < C ^ ~2.7, i.e. the DM 
profile in galaxy clusters either follows the NFW profile or 
is slightly steeper at large radii; 

• for systems at redshift z = 0, with NFW DM, F = 6/5 
and Too = 0, we also study the dependence of the theo- 
retically derived M-T relation on concentration and find a 
very weak dependence, unless the temperature is calculated 
without excising the central region; 

• for systems with NFW DM, concentration given by 
DOS, F = 6/5 and Too = 0, we study the dependence of the 
theoretically derived M-T relation on redshift, and find that 
our model scales in a self-similar behaviour. However, given 
the freedom in our model, the departure from self-similarity 
seen in some clusters can be used to constrain the evolution 
of the polytropic index with redshift (see e.g. Equations 1261 
and|30j; 

• for systems at redshift 2 = 0, with NFW DM and Too = 
0, we study the dependence of the theoretically derived L- 
T relation on the polytropic index, on concentration and 
on the stellar mass fraction. We find that a concentration 
strongly dependent on mass is ruled out by comparison to 
recent X-ray observations. 

The power of this model lies in its simplicity and its 
flexibility. In future work, more sophisticated models can be 



built upon this framework by relaxing some of the assump- 
tion made here, e.g. spherical symmetry, and by including 
non-thermal contributions to HE; a radially-dependent poly- 
tropic index; more general DM profiles; gas self-gravity; and 
the effect of the stellar gravitational potential. We also note 
that we studied narrow slices of the five-dimensional param- 
eter space (F, To, po, C a-nd Cvir). Further studies can easily 
exploit the model to explore the full parameter space. 

Although we have used the polytropic relation to de- 
scribe the gas, its origin is still unclear, although it is likely 
tied to the formation of the halo. In this paper, we have 
proposed a relation between the polytropic index, mass and 
redshift, thus providing useful information and a first step 
to solve this puzzle. 
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Figure Al. DM radial profiles for a system of virial mass Mvir = 
IO^'^Mq at redshift z = 0, with concentration given by D08, for 
three cases of the BIO profile: fS = 3/2 {black, dotted curves), f} = 
2 (the NFW case; black, solid curves) and /3 = 5/2 (black, dashed 
curves). The black, vertical, dot-dashed lines denote the virial 
radius, whereas the black, vertical, dashed lines denote the DM 
scale radius. Top panel: DM density (thin curves) and average 
density (thick curves) radial profiles. The horizontal lines denote 
the A vir /dm Pc (black, solid), 200/dmPc (blue, dotted), 500/dmPc 
(red, dashed) and 2500/dmPc (black, dot-dashed) density values, 
respectively. Bottom panel: DM enclosed mass radial profiles. 
The black, horizontal, solid line denotes the DM virial mass. 



APPENDIX A: 
RADII 



OVERDENSITY MASSES AND 



In this section, we describe the difference between an NFW 
profile and a BIO profile with /3 7^ 2. In Figure KT] we show 
the DM density pdm, average density p^j^j and enclosed mass 
radial profiles for a system of virial mass Mvir = IO^^Mq at 
redshift z = 0, with concentration given by DOS, for three 
cases of the BIO profile. By definition, the virial radius r^ii is 
the same, regardless of the choice of /3. If we assume that the 
concentration Cvir does not change for any 1 < /? < 3, the 
scale radius Vs is the same for all profiles. What changes in- 
stead is the characteristic density S^, as described in Section 

m 

The DM average density and enclosed mass are higher 
(lower) than the NFW case in the /3 > 2 (/3 < 2) case for 
r < rvir, and lower (higher) than the NFW case in the /? > 2 
(P < 2) case for r > rvir. This results in different overdensity 
masses and radii; for /? > 2 (/? < 2), Ma^, and are greater 
(smaller) than in the NFW case, for any Ac > Avir (which 
is the case at any redshift for Ac > 200). 
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Figure A2. DM overdcnsity masses (top panel) and radii (bot- 
tom panel), for systems at redshift z = 0, with concentration 
given by DOS, for three cases of the BIO profile (/3 = 3/2 - dotted 
lines; = 2 — the NFW case, solid lines; and /3 = 5/2 - dashed 
lines) and three values of Ac (Ac = 200 - blue lines; Ac = 500 
- red lines; and Ac = 2500 - black lines). The black, dot-dashed 
line in the bottom panel is the virial radius. 



For a system of known Mvir and known concentration, 
assuming that the baryonic mass is negligible, we can calcu- 
late the overdensity mass Ma^ and the overdensity radius 
tAc, for any given value of Ac and p. By combining the def- 
initions of virial mass and overdensity mass with those of 
enclosed mass, we obtain rA^ by numerically solving 



d Ac 



+ [(1 - 0){rAjrs) - + taJts)'-^ 



r^i,Avi, 1 + [(1 - 

in the /3 7^ 2 case and 



)Cvir - 1] (l + Cvi,)'"" 



rj^Ac _ In (1 + TAjr,) - {rAjr,)/{l + r^^/rj 



In (1 -I- Cvir) - Cvir/(1 + Cvir) 

2 case. We subsequently solve 



7-A Ac 



(Al) 



(A2) 



(A3) 



T-3 A ■ 
' vir"vir 

in the j3 

Ma, 

Mvir ^ J-virAvir 

to obtain Ma,. 

In Figure IX2I we show the overdensity masses and radii 
as a function of virial mass, for four typical values of Ac 
(including Avir) and three values of /3, when using the mass- 
concentration relation by DOS at redshift 2 = (when 
Avir ^ 95). 
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3/2 


Ac = Avir 


Ac = 200 
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Ml 
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308514 
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-0 591 205 
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5/2 


A A 


— zUU 


ZAc — OUU 


ZAc — ^OUU 


Ml 







-0.017403 


-0.037090 


-0.058379 


Ma 
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0.871441 


0.704620 


0.408453 


Ri 




0.333333 


0.330437 


0.325672 


0.312246 


R2 




0.080661 


-0.046555 


-0.210196 


-0.523894 



Table Al. Overdensity fitting parameters. We list the fitting 
parameters for Equations l)A4[l and l|A5|l . for systems at redshift 
z = 0, with concentration given by D08, for four typical values of 
Ac - 200, 500, 2500 and Avir(z = 0) ~ 95 - and for three values 
of /3 - 3/2, 2 (NFW) and 5/2 - in the BIO profile, when ta^ is 
given in Mpc and A/vir is given in units of 10"'^*Mq. 



Given the weak dependence of concentration with mass 
(DOS), we can fit simpl(|3 relations between Ma,, ^Ac and 

Mvir: 



Ma, 

Mvir 

and 



= Mllogio(Mvir) + M2 



logio(rAj = R\ logio(A/vir) + i?2. 



(A4) 



(A5) 



where the parameters Mi , Ma , R\ and R2 depend on Ac and 
/3, as shown in Table IXT1 for the z = case, when ta, is given 
in Mpc and Afvir is given in units of IO^^Mq. i?i_vir = 1/3 
simply from the definition of virial radius. 



APPENDIX B: THE COOLING FUNCTION 

Cooling, defined as the total energy lost per unit volume per 
unit time, can be usually parametrized in the optically-thin, 
fully-ionized, low-density limit as Anct = nirieAN, where the 
cooling function An is defined below. We assume, following 
Tozzi & Norman (2001: TNOl), n-.n^ = 0.704 (p/ryip)^ for a 
gas with mean metallicity [Fe/H] = —0.5. 

Cooling time, defined as the total energy of the gas di- 
vided by the energy loss, is then 



''"cool 



(ric +ni)fcBT 



(Bl) 



(7 - l)ncniAN n* (7 - 1)An ' 

where n« = noni/(nc -f ni) and we have assumed that elec- 
tron and ion temperatures are equal to each other. 

If we use mass-concentration relations with a stronger mass 
dependence (e.g. OlO and Oil), the fitting relations are not this 
simple. 
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[dimensioiiless] 


Ca 







1.449 




f \ 1 TYir>TlCir\Tl It^GG 


^5 


Q 








[kcV^"*"] 


^0 







6.303 




r\ 1 TYi r>n Qi mi lr>c;c; 

11 1*^11 It^OO 


C7 







5.0 X 10~^ 




[keV~^] 


c; 
^1 


8.8 X 10" 


13 


2.256 X 10" 


14 


[erg cm^ K ] 




1.7 X 10" 


27 


1.540 X 10" 


27 


[erg cm^ ] 


r' 


6.3 X 10" 


24 




24 


t;l t L,lll 


a' 


-1.7 




-1.443 




[dimensioiiless] 




5.0 X 10" 


1 


5.0 X 10"^ 




[dirQensionless] 









1.449 




[dimensioiiless] 









4.551 X 10" 


6 


[K-i] 









6.303 




[dimensioiiless] 









4.309 X 10~ 


10 


[K-i] 



Table Bl. Cooling function fitting parameters. We list the fitting 
parameters for Equation (BS), for a mean metallicity of [Fe/H] = 
—0.5, for both the TNOl and this paper cases. When substituting 
k^T (in keV) with T (in K), the primed parameters in the table 
should be used instead. 



Figure Bl. Cooling function (top panel) and normalised cool- 
ing time (bottom panel) for an optically-thin, fully-ionized, low- 
density gas of mean metallicity [Fe/H] = —0.5 in coUisional ion- 
ization equilibrium. Black, solid lines: data from Sutherland & 
Dopita (1993). Blue, dotted lines: fitting function to the data by 
Sharma et al. (2010). Red, dashed lines: our fitting function to the 
data. The two black, vertical, dotted lines denote the temperature 
regimes of Equation IIB3I1 . 



keV"-' ^^ Bl = 7.027 x 10"''' erg cm'^ s"' keV""=, A2 = 6 
and B2 = 0.6. 

We adopt instead a new fitting function, wtiicii iias tiie 
advantage of being a single formula like the one by TNOl, for 
easier computational and/or analytical use, and at the same 
time has the advantage of taking into account the decrease 
of An at low temperatures (for log\Q(r[K]) > 4.45) and the 
relativistic correction for the Bremsstrahlung emission: 



Under the assumption of coUisional ionization equilib- 
rium, we adopt the cooling function by Sutherland & Dopita 
(1993). TNOl first provided the fitting function, valid only 
for ksT > 0.02 keV, 



An = CiiksT)" + C2{kBTf + C3, 



(B2) 



where the constants Ci, C2 and C3 depend on the mean 
metallicity of the gas, a — —1.7 and /3 = 0.5. For a gas with 
mean metallicity [Fe/H] — —0.5, the constants are Ci = 
8.6 X lO"^'^ erg cm^ s" ' 



keV"°, C2 = 5.8 X 10"^'' erg cm^ 



"' keV"'', and C3 = 6.3 x 10"^" erg cm^ s" 



Sharma, Parrish, & Quataert (2010) recently gener- 
alised the fitting function of TNOl in the [Fe/H] = -0.5 
case, by considering three separate temperature regimes: 



Ci(fcBr)°exp 



C4 



1 + CsfcsT 



+ 



(B4) 



+ C2{kBTf{l+CTkBT)+C3. 



In Figure IBll we show the cooling function and nor- 
malised cooling time, n,Tcooi, as a function of temperature, 
for the [Fe/H] = —0.5 case. For comparison, we show the 
original data from Sutherland & Dopita (1993), the fitting 
function by Sharma et al. (2010) and the new fitting function 
introduced in this paper. 

In Table [BTI we list the fitting parameters for Equation 
()B4|I . for a mean metallicity of [Fe/H] — —0.5, for both the 
TNOl and this paper cases. 



A^ 



'Ai(fcBr)^^ 

BiikBTf^ 



if r < Ti, 
ifT2>T>Ti, (B3) 



CiiksT)" + C2ikBTf + C3 ilT>T2 



where /cbTi = 0.0017235 keV, kBT2 = 0.02 keV, the con- 
stants Ci, C2, C3, a and /3 are the same as in TNOl and 
the other constants are Ai = 5.891 x 10~® erg cm^ s~' 
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